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In this paper, we investigate the stationarity of stochastic processes in the fractional Fourier domains. 
We study the stationarity of a stochastic process after performing fractional Fourier transform (FRET), 
and discrete fractional Fourier transform (DFRT) on both continuous and discrete stochastic processes, 
respectively. Also we investigate the stationarity of the fractional Fourier series (FRFS) coefficients of a 
continuous time stochastic process, and the stationarity of the discrete time fractional Fourier transform 
(DTFRFT) of a discrete time stochastic process. Closed formulas of the input process autocorrela- 
£> . tion function and pseudo-autocorrelation function after performing the fractional Fourier transform are 

derived given that the input is a stationary stochastic process. We derive a formula for the output 

' th 
\q ' autocorrelation as a function of the a power spectral density of the input stochastic process, also we 

\Q \ derived a formula for the input fractional power spectral density as a function of the fractional Fourier 

O 

. transform of the output process autocorrelation function. We proved that, the input stochastic process 

must be zero mean to satisfy a necessary but not a sufficient condition of stationarity in the fractional 
domains. Closed formulas of the resultant statistics are also shown. It is shown that, in case of real 
input process, the output process is stationary if and only if the input process is white. On the other 
hand, if the input process is a complex process, it should be proper white process to obtain a stationary 
output process. 

Index Terms 

Fractional Fourier transform (FRET), discrete FRFT (DFRFT), additive white Gaussian noise(AWGN), 
fractional power spectral density 
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I. Introduction 

The fractional Fourier transform (FRFT) has been used in most recent applications and 
frequently used as a tool in signal processing and analysis. It has been discussed in many papers 
and has been proved to be very useful in solving some problems in quantum physics, optics, and 
signal processing [JTJ— [JTTJ. In 021, the authors studied the stationary of continuous signals in 
the fractional domain. The relationship among the FRFT, the linear canonical transform (LCT), 
and the stationary and nonstationary random processes is derived. Surprisingly, they found many 
interesting properties. For instance, if we perform the FRFT for a stationary process, although the 
result is no longer stationary, the amplitude of the autocorrelation function is still independent 
of time. For the FRFT of a stationary process, the ambiguity function (AF) is a tilted line and 
the Wigner distribution function (WDF) is invariant along a certain direction. In the same paper, 
they used the notion of the fractional stationary random process and found that a nonstationary 
random process can be expressed by a summation of fractional stationary random processes. 

The purpose of this paper is to study stationarity of stochastic processes in the fractional 
Fourier domains. We study stationarity of stochastic processes after performing the FRFT and 
DFRFT. We also investigate the stationarity of the fractional Fourier series (FRFS) coefficients 
of a continuous time stochastic processes and the stationarity of the discrete time fractional 
Fourier transform (DTFRFT) of a discrete time stochastic process. The case of complex Gaussian 
processes and circular symmetric Gaussian processes were investigated. 

The effect of performing Fourier transform on stationary proper process was studied ffT3l . It 
has been shown that the discrete Fourier transform (DFT) of the input sequence is proper if and 
only if the input sequence is proper. The term proper refers to a process x(t) with C(ti,t 2 ) = 0, 
i.e., complex random variables and processes with a vanishing pseudo-covariance are called 
proper, where 

C(h,t 2 ) = EUxih) - £{z(ti)})(z(t 3 ) - E{x(t 2 )})\ (1) 

is the relation function (pseudo) which is necessary for complete description of second-order 
statistic Ifl4l . 

If the discrete time stochastic process z[n] is a wide sense stationary (WSS) process with power 
spectral density Z(u), then its DFT, Z(ui), is a nonstationary white noise with autocovariance 
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is given by [15, p. 519]: 

El Z(u)Z*(v) )■ = 2irS z (u)5(u - v), -n < u, v < tt 



Z(u)= z[n]e- inuJ (2) 

n=— oo 

z[n} = — f Z(u)e inuJ duj 

where u = 2irf, f is the frequency domain index and S z (u) is the power spectrum of z[n]. 

In lfT3l . it has been shown that the properness is preserved under affine transformations and 
also the complex-multivariate Gaussian density assumes a natural form only for proper random 
variables. It was proved that circular stationarity of a proper complex time-domain sequence 
z[0], z[l], . . . , z[N— 1] corresponds to uncorrelatedness of the components of its DFT (frequency- 
domain sequence Z[0], Z[l], . . . , Z[N— 1]), i.e., the autocorrelation function of the DFT sequence 
is given by: 

E^Z(k)Z*(£)^ = VNS z (k)5(k - £) (3) 

where k, i E {0, 1, . . . , N — 1}, S z (k) is the DFT of the input sequence autocorrelation function, 
and N is the input sequence length. 

The main contributions of this paper can be summarized as follows. 

. We investigate the stationarity of a stochastic process after performing FRFT, DFRFT, FRFS 
and DTFRFT. 

• We study general features of the output signal, i.e., mean, autocorrelation, and pseudo- 
autocorrelation for the discrete and continuous cases. We provide closed formulas of the 
output features. 

• We investigate stationarity of real and complex stochastic processes. In addition, we provide 
necessary conditions on the input process for stationarity of the output process. 

• The relationship among the autocorrelation, pseudo-autocorrelation and the power spectral 
density of a process is derived. 

This paper is structured as follows. In Section lU we study stationarity in continuous FRFT. In 
Section Unl we investigate the stationarity of the FRFS coefficients of a continuous time stochastic 
process and the DTFRFT of a discrete time stochastic process. The stationarity, statistics and 
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probability density function after performing DFRFT is discussed in Section [IV] and finally, we 
conclude the paper in |V] 

II. Continuous Fractional Fourier Transform 

The continuous FRFT (Fig. Q]) is an integral transformation with a — a- where a represents 
the fractional order. The FRFT of a function x(t) is given by: 



/oo 
z{t)K a (t, 
-oo 



u)dt (4) 



where Ta^u denotes a transformation from t to u with angle a, defined as the rotation angle 
with respect to t axis, and K a (t,u) is the transformation kernel defined as following: 



K a (t,u) 



2 

1— i cot a gi cot a 



2- 



-itucsca+i^- cot a :f _L _„ 

5(w — i) if a = 27rp 

5(u + t) if a + 7r =2np 

where p = — oo, . . . , — 1, 0, 1, . . . , oo. The fractional Fourier transform with order a is a rotation 
of the time-frequency plane with respect to time axis by an angle a = |a, the Fourier transform is 
a rotation by an angle a — | Q. The FRFT kernel can be represented by its eigendecomposition 
functions |[T6ll 

oo 

K a (t,u) = Y,e- ian H n (t)H n (u) 

n=0 

where H n (.) is the normalized Hermite function with unitary variance. The inverse formula of 
the fractional Fourier transform is obtained through replacing a by —a as follows: 

/oo 
Z a (u)K- a (t,u)du 
-oo 

In this paper, we assumed a takes any value between — | and |, excluding, the special cases, 
where a = 0, ±|. We assumed that the FRFT is a system with two terminals the input terminal 
and the output terminal. In Eqn.(|4]), the input is z{t) and the output is Z a {u). Assume that z(t) 
is a stochastic process with autocorrelation function given by [TT3l . 031, tfTTl : 

ft,(*l,*2)=£{*(*l)**(*2)} (6) 

and a pseudo-autocorrelation function: 
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K z {t 2 M) = E{z{h)z{t 2 )} (7) 

where E{.} denotes the expectation, and z(t\), z(t 2 ) are two samples of the original stochastic 
process z(t), taken at time instants ti and t 2 , respectively. The expected value of the output 
process, after performing FRFT, is given by: 



fi a (u) = E{Z a (u)} = / E{z(t)}K a (t,u)dt 

J — oo 

/I a (u) = Fa,t^u{lJ>{t)} (8) 

where p(t) = E{z(t)}, and, fx a (u) is the FRFT of the input process mean fi(t). If the process 
is wide sense stationary, fi(t) is independent of time, = ji, then the expected value of the 
output process is: 

H a (u) = [I x 7 a ^ u {\) (9) 



^„{i}=v^ ± ^ e " i * , "° (10) 

It can be shown that from (|8l.((9l).(fT2l) the expectation of the output of FRFT is non-stationary, 
i.e., output mean depends on u, even if the input process is wide sense stationary. If the input 
process is zero mean, ji — 0, the expected value of the output process is given by: 

M«) = o. (11) 



Proposition 1: Given a wide sense stationary process z(t) with mean, fx ^ 0, the output of 
the fractional Fourier transform is a non- stationary process with mean fx a (u) 



„„(„) = ^i±l^e-*— (12) 

Since a time variant mean is a necessary condition of stationarity, therefore, we obtain the 
following. 

Proposition 2: A zero mean input process is a necessary but not a sufficient condition for 
the stationarity in the fractional domain. 

Proof: From Proposition \T\ the expected value in the fractional domain fi a (u) is always 
a function of (varying with) the fractional index u except when /i = 0, therefore, the output 
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process violates a necessary condition of stationarity which is the independency of the output 
process on u. ■ 
It should be mentioned that complex random processes with nonzero mean are usually not of 
interest, since a complex envelope with nonzero-mean corresponds to a non-stationary band- 
pass process, i.e., the expected value of the bandpass process of non-zero mean with mean 
baseband signal is E{^t{z{t)e l2 ' nfct }} = §l{E{z(t)}e l7T - fct } which is a time dependent process, 
since E{z(t)} ^ [13]. Let us use equation ([4]) to derive the autocorrelation function of the 
output of the FRFT as a response to an input process z(t) described above ©. Assuming 
J^a^uizit)} = Z a (u), the output autocorrelation function is given by: 

K a ( Ul ,u 2 ) = E{Z a ( Ul )Z* a (u 2 )} (13) 

/OO POO 
/ E{z(t)z*(s)} 
-oo J —oo 

K a (t, u^K^s, u 2 )dsdt 

TZ(s, t)K a (t, ui)K^(s, u 2 )dsdt 



oo roo 



oo J — oo 



oo 



1Z(s,t)K a (t,u 1 )dt 



K* a (s,u 2 )ds. 



/oo 
(3 a (ui,s)K_ a (s,u 2 )ds. (14) 
-oo 

where /3 a (ui, s) = lZ(s, t)K a (t, u\)dt and it represents a rotation of the input autocorrelation 
function by angle a. The equations above can be interpreted as follows. The resultant autocorrela- 
tion function is a rotation of the input autocorrelation function with respect to t by angle a given 
certain s, then rotation of the resultant of the first rotation with respect to s by angle —a given a 
fixed u\. The same argument can be shown for the output pseudo-autocorrelation function. The 
relationship among the autocorrelation, pseudo-autocorrelation and the power spectral density of 
a process is derived in Appendix B. 
Since the Kernel is unitary, i.e., 

K a (s,ui) K^(s,u 2 )ds = 5(u 2 — u-\), it can be shown that the output autocorrelation 
function as a response to an input autocorrelation function 1Z(t 2 — ti) = ^£-5{t 2 — ti) is given 
by: 

K a {u 2 - ui) = -£-8{u 2 (15) 
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z(t) 



FRFT 

a 



Fig. 1. FRFT block 

Furthermore, the output pseudo-autocorrelation function as a response to a proper and zero 
mean process is given by: 



Therefore, properness is preserved under the fractional Fourier transform. The proofs of (|T"5) 
and (fT6l) are existed in Appendix C. 

Proposition 3: Given a real wide sense stationary process z(t) with mean /i and autocorrela- 
tion 1Z(t) the output process of the FRFT Z a {u) is a non- stationary process. 

Proof: Since lZ a (ui,it2) is a function of domain indexes u\ and u 2 , therefore, the output 
process Z a (u) is a non- stationary process and stationarity can be obtained if and only if (iff) 
the input process is a white process (See Appendix A and C). ■ 

Proposition 4: Given a complex wide sense stationary process z(t) with mean /i, autocor- 
relation 7Z(t), and pseudo-autocorrelation 7Z(t) the output process of the FRFT Z a {u) is a 
non-stationary process and stationarity can be obtained iff the input process is a proper white 
process. 

Proof: Since lZ a {ui,U2) is a function of domain indexes u\ and u 2 , therefore, the output 
process Z a {u) is a non- stationary process and stationarity can be obtained iff the input process 
is a proper white process (See Appendix A and C). ■ 

III. Discrete time fractional Fourier Transform and Fractional Fourier 

Series Expansion 

In this section, we study stationarity of the discrete time fractional Fourier transform (DT- 
FRFT) and fractional Fourier series (FRFS). Both DTFRFT and FRFS were studied in lfl8l . 
Aperiodic signal z(t) on a finite interval t E [— T/2, T/2] can be expanded by its FRFS as 



TL a {u 2 - Mi) = 0. 



(16) 
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following: 

z(t)=J2 C a,n(j>a,n( t ) ( 17 ) 
n 

where n = — oo, . . . , — 1, 0, 1, . . . , oo and <f) a n is orthonormal basis and is given by: 

= (18) 

'Tcsca 



2vr 

where t Q is called the central frequency in the fractional Fourier domains. However, t a must 
equal to 2ns ™ a , to guarantee that (f) a>n (t) satisfy the orthogonal condition. 

The FRFS expansion coefficients are computed by the inner product of the signal and chirp 
basis signals 4> a ,n{t)- In Ifl8lh it was proven that the FRFS coefficients of an aperiodic signal 
z(t) on a finite interval t G [— T/2, T/2] can be obtained from the sampled values of FRFT and 
it is given by: 

C v =^p°z a (^pL) (19) 

where Z a (.) is the FRFT of z(t), and T is the interval width. 
The expected value of C a n is given by: 



nf/1 , / 2ix sin a _ , _ . 27rsina. 1 

E{C ain } = \j^—E{Z a {n^—)}. (20) 

Same argument as in Section HIT] about the dependency of the expected value of C a>n on the 
fractional index n can be established here. Therefore the coefficients are not stationary in general, 
only if the input stochastic process is zero mean, the output mean will be independent of n with 
zero mean. The autocorrelation function of the FRFS coefficient C Qjn is given by: 

E{C^C- J)} = ^ E{ z a{ n^)Z- aV ^)} (21) 

where n, i = — oo, . . . , —1, 0, 1, . . . , oo. The autocorrelation function of two coefficients of the 
FRFS can be stated as: 

E{C^ /)} = ^p Ka(n ^, (22) 

Similarly, we can show that the pseudo-autocorrelation function is given by: 

i^c^ ^ \ i 27rsina - . 27rsina „27rsina. 

E{C a , n C ai£ )} = — Tl a {n — ,£ —— ). (23) 
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Sampling a continuous time bandlimited signal z(t) by sampling time T s , we obtain z[n}. The 
DTFRFT D a [k] of z[n] is related to FRFS coefficients as following lfT8l : 

D a [k] = C^ /2+a>k . (24) 

where k = — oo, . . . , — 1, 0, 1, . . . , oo. 

Proposition 5: Given a finite support wide sense stationary process z(t) with < t < T, 
mean \i, autocorrelation function TZ(r) and pseudo-autocorrelation function 7t(r), the FRFS and 
DTFRFT coefficients of z{t) are non-stationary processes. 

Proof: Since the FRFS and DTFRFT coefficients of z(t) are functions of TZ a (ui,u 2 ) which 
is a function of domain indexes U\ and u 2 , therefore, the coefficients are non-stationary processes 
and stationarity is obtained iff the input process the input process is white in case of real process 
and proper white process if the input process is complex (See Propositions |3] and H]). ■ 

IV. Discrete fractional Fourier Transform 

In this section, we are studying the effect of performing the DFRFT (Fig. [2]) on a discrete 
stochastic vectors. The DFRFT is an affine unitary transformation. It is the eigendecomposition 
of the ordinary DFT matrix, i.e., T a = QA a Q H , where, Q is the eigenvectors matrix of the 
DFT matrix, and A is the eigenvalues matrix. Assuming that the sampling process produced A" 
samples, i.e., time-bandwidth product equals N. 

Define a stochastic process {Zi} where i = 1,2, Let us rearrange the process samples 

as following, Z = Z 2 , . . . , Z N ] T , where, Z is N x 1 complex column vector, and Z k = 
Xk + iYk is a stochastic random variable where XkXk G M, 2 = with mean [i = E{Z} = 

/i(2), . . . , /i(AQ] T , covariance matrix C z = E{(Z — /S)(Z — jj) H }, and pseudo-covariance 
(relation matrix) P = E{(Z — ji){Z — jj) T }, which is necessary for complete description of 
second-order statistic [[141 [[191 where T and H denote transpose and hermitian (transpose and 
conjugate), respectively, performing discrete FRFT on process vector Z the output process Z a 
is given by: 

Z a = 7 a Z 

The output Z a is a random vector with mean \i a , pseudo-covariance matrix P a and covariance 
matrix C a 

^ a = E[Z a ]=E{T a Z} = Ty. (25) 
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Equation (1231) can been seen as the discrete FRFT of /i. The discrete FRFT matrix maps N 
points from time domain to N points in fractional domain, each sample will be mapped with 
different mean based on its location in time and fractional domain, i.e., the expected value of 
the nth output equals the inner product of the nth row of the discrete FRFT matrix and n H . The 
covariance matrix C a is given by: 

C a = E^(Z a -u a )(Z a -fi a ) H ^ 

= J*e{{Z-h){Z-h) b }t- 
Covariance matrix of Z is C z hence, 

C a = T a C z T~ a . (26) 
Similarly, the output pseudo-covariance matrix will take the form: 

P a = F a P[f a f. (27) 

Since T a is affine linear unitary transformation, hence, when the input is Z ~ CJ\f(u,,C z ,P), 
the output of the discrete FRFT is Z a ~ CJ\f(fx a , C a , P a ), where CJ\f(.) denotes complex normal 
distribution. For iV-dimensional complex Gaussian random vector Z has real Z^ and imaginary 
Zq. components which form a 2iV-dimensional real Gaussian random vector and Z is proper 
i.e. its pseudo-covariance matrix vanishes lfl3l . the probability density function (pdf) of Z a ~ 
CN{na, C a ) is given by [HI, [[H 



t (7\— _ r -{Zg-v. a ) H C a 1 {Z a - l i a ) 

fzA a) ~ n N det(C a ) 1 } 

If Z is proper with zero mean, the covariance matrix fully specifies the first- and second-order 
statistics of a circular symmetric random vector [H9l . If Z ~ CJ\f(0, <j 2 Inxn) which is a circular 
symmetric Gaussian random variables, the output of the discrete FRFT is Z a ~ CJ\f(0, (t 2 Inxn)- 



V. CONCLUSION 

In this paper, we have discussed the stationarity of a stochastic process in fractional domains. 
We have found that the zero mean input process is a necessary but not sufficient condition for 
stationarity in the fractional domains, and that the fractional domains are non- stationary. It has 
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T >Z a (3) 

J ^(tf) 



Fig. 2. DFRFT block 



been shown that, in case of real input process, the output process is stationary if and only if the 
input process is white. On the other hand, if the input process is a complex process, it should be 
proper white process to obtain a stationary output process. The statistics of the output stochastic 
process were proved and discussed. The results are a backbone in systems based on fractional 
Fourier transform such as (OFDM) systems [0, sonar, radar systems, moving target detection in 
airborne SAR [ETI . chirp detection [|22l . [|23l . speech processing JT), Optical image encryption 
0, beamformers [|24], and many other applications in signal processing and optics [|4]|. 

Appendix A 

In this Appendix, we derive closed formulas of the output autocorrelation function and pseudo- 
autocorrelation function corresponding to a stationary input process. Assume that the input 
process is stationary process with autocorrelation function 1Z(s — t) where s,t are two instants 
of the input process z, the output autocorrelation function is given by: 

TZ a (u u u 2 ) = E{Z a (ui)Z*(u 2 )} 

K(t - s)K a (t,u 1 )dt 



K* a (s,u 2 )ds 

S z>a (ui -cos(a)s)e misin(a)s 

) 

e islnia)cos{a) ^K* a {s,u 2 )ds 

where S z , a {ui) = ^ 7 Q,t->u 1 {'^(t)} and it represents the FRFT of the input autocorrelation function 
with respect to t (a th fractional power spectral density), the output is a transformation from t to 
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u\ with order a. Using properties of the fractional Fourier transform [|5], the autocorrelation as 
in lfT2l is given by: 



TZ a (u 1 ,u 2 ) = sec(a)^(sec(a)(n 1 - U2 )y^ 2 - Ul )tM*) (2 9) 

where a ^ , t = 1, 3, 5, . . . , oo. Similarly, we can find the output pseudo-autocorrelation 
function as a response of the stationary input process, according to ITPJI the wide sense stationary 
process has pseudo-autocorrelation function depends on time shift, i.e., independent of time 
Tt(r) = 7Z(t — s). With some change of variable we get: 

n a ( Ul , u 2 ) = sec(a)T a , y ^ U2 |g q (u x -y cos a ^ {ui '" 2)y | (30) 

where 

c = sin a cos a + cot a 



T a (ui,u 2 ) = (u\ sin a — M2Csca)W h -^csca. 

cot a 

performing some properties of the FRFT (i.e., shifting, scaling and exponential multiplication 
properties), and define (3 = arctan(c 2 tan a), we get the following: 



K a {ui, u 2 ) = \ — : — — e 

cr — % cot a 



1-icota ^(^Acos^) 

v 1 \ cos z a / 



A sm 3 

G 2 b((u 2 - r^(Mi,M 2 ) sina) — : h «i cos/3) 

v c sm a ' 

g— iu\ (u2— («i ,U2) sin a) sin /3^j sin fi cos (31) 
g— « sui(o) cos(a) 2 
e «r a (ui,U2)u2 cos(a) 



where T a . a ^ u A Gp(ui) \ = G 2f} {u 2 ) = T 2 p,t-+u % \ Tl{t) >■ 



Appendix B 

In this Appendix, we prove the relationship between the output process autocorrelation function 
and the fractional power spectral density. The power spectral density of a stationary process is 
the Fourier transform of the autocorrelation function, a th fractional power spectral density is the 
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fractional Fourier transform of the autocorrelation function with order a. The autocorrelation 
function of the output of the fractional Fourier transform is given by: 

TZ a (u 1 ,u 2 ) = E^ Z a (ui)Z^(u 2 ) > 

J (32) 



lZ(t—s)K a (t, ui)K^(s, u 2 )dsdt. 
Let us substitute by K a (t,ui) and define r = t — s the resultant is: 



00^^30 



1 — i cot a 



OO poo 2 



00 J —00 

2 

— i csc(a)rui — t csc(a)sui cot(«) 
6 1 



2 9 
i^cotia) i^- cot(a) 



e 2 v ; e 2 

K*(s, u 2 )dsdt. 
e 4^K*(s,u 2 )ds 

define S z>a (ui) as the a th fractional power spectral density of the input process z, i.e., 

S z , a (Ul) = J = a,r^u 1 {^i T )} 

POO 2 

K a { Ul ,u 2 )= / 5 Zia (n 1 -cos(a)s)e- icot2(a)sin(a)cos(a) Ve iUlCOt ^ cos ( Q ) s 



(33) 



e -i csc(a) S «i cot(a) 



.2 



^(tirt) = ^ Q ,™<! ^, a (ui - cos(a)s)e tcos ^ sia ^e- is ' m ^ s }. (34) 



With the same procedure as in [12], we can express 1Z a (ui,u 2 ) as a function of the fractional 
Fourier transform of the fractional power spectral density, define J r _ a ,j,-> lt2 {5 , 2)Q ,(y)} = lZ(u 2 ) it 
can be shown that: 

S z , a (ui ~ cos(a)s) = Je-^W^^e*^)" 1 - (35) 

by making the change of variable, u = U\ — cos as, the fractional power spectral density can be 
expressed as a function of the output autocorrelation function as follows: 

S 2 , Q (u;)=J^|ft Q (a;+cos^ 

where 1Z a (.) is given in Appendix A. 
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Appendix C 

In this Appendix we will prove the autocorrelation function of the output process, when 
the input is a white stochastic process, i.e., additive white Gaussian noise (AWGN), with 
autocorrelation function TZ(r) = ^6(r), where r = s — t is a time shift, N a /2 is the power 
spectral density of the process: 

7fc a («2-tti) = E{Z a {u 1 )Z* a {u 2 )} 



OO /"00 



oo J — OO 



E{z(t)z*(s)}K a (t, Ul ) 
K*(s, u 2 )dsdt 

N f°° f°° 

— / K a (s,u 1 )K*(s,u 2 ) / 5(s — t)dtds 

" J —OO J —OO 

S v ' 

=lVs£[— 00,00] 

K a (s,ui)K* a (s,u 2 )ds 



No 
2 



oc 



Kernel is unitary, hence, 

/OO 
K a (s, u 2 )ds = 5(u 2 - m). 

-00 

Finally, the resultant autocorrelation function 



N 

TZ a (u 2 - ui) = — -8(u 2 - ui). (36) 



This proof can be obtained through substituting in Eqn. (1291 ) by the value of the input process 
autocorrelation function and notifying that sec(a)<5(sec a(u\ — u 2 )) = 5{u\ — u 2 ). 

Similarly, the output pseudo-autocovariance function as a response to a proper stationary input 
process is given by: 



00 roo 



00 J —00 



C a {u u u 2 ) = E\ (Z a (ui) - ii<x){ux){Z a {u 2 ) - pL a (u 2 )) 

E{(z(t) - ii){z{s) - [i)}K a (t, ui)K a {s, u 2 )dsdt 
= 0. (37) 

Hence, the properness is preserved under the fractional Fourier transformation, for the WSS, 
proper, and zero mean process described in (|36l . the pseudo-covariance converges to the pseudo- 
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autocorrelation of the process, and, both converges to zero: 

C a (u 2 , Ml) = TZa(u 2 , «l) = 0. (38) 



References 

[1] D. Ma, X. Xie, and J. Kuang, "A novel algorithm of seeking FrFT order for speech processing." Prague, Czech Republic: 

IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), May 2011, pp. 3832-3835. 
[2] R. Tao, J. Lang, and Y. Wang, "Optical image encryption based on the multiple-parameter fractional Fourier transform," 

Optics letters, vol. 33, no. 6, pp. 581-583, 2008. 
[3] Q. Yang, R. Tao, Y. Wang, and E. Chen, "MIMO-OFDM system based on fractional Fourier transform and selecting 

algorithm for optimal order," Science in China Series F: Information Sciences, vol. 51, no. 9, pp. 1360-1371, 2008. 
[4] H. Ozaktas, M. Kutay, and Z. Zalevsky, The fractional Fourier transform with applications in optics and signal processing. 

Wiley New York, 2001. 

[5] L. Almeida, "The fractional Fourier transform and time-frequency representations," IEEE Transactions on Signal 

Processing, vol. 42, no. 11, pp. 3084-3091, 1994. 
[6] A. Zayed, "A convolution and product theorem for the fractional Fourier transform," IEEE Signal Processing Letters, 

vol. 5, no. 4, pp. 101-103, 1998. 
[7] H. Ozaktas, O. Arikan, M. Kutay, and G. Bozdagt, "Digital computation of the fractional Fourier transform," IEEE 

Transactions on Signal Processing, vol. 44, no. 9, pp. 2141-2150, 1996. 
[8] T. Alieva, V. Lopez, F. Agullo-Lopez, and L. Almeida, "The angular Fourier transform in optical propagation problems," 

Journal of Modern Optics, vol. 41, pp. 1037-1044, 1994. 
[9] D. Mendlovic and H. Ozaktas, "fractional Fourier transforms and their optical implementation: I," Journal of the Optical 
Society of America A, vol. 10, no. 9, pp. 1875-1881, 1993. 
[10] A. Lohmann, "Image rotation, wigner rotation, and the fractional Fourier transform," Journal of the Optical Society of 

America A, vol. 10, no. 10, pp. 2181-2186, 1993. 
[11] V. Namias, "The fractional order Fourier and its application to quantum mechanics," Journal of the Institute of Mathematics, 

vol. 25, no. 10, pp. 241-265, 1980. 
[12] S. Pei and J. Ding, "Fractional fourier transform, wigner distribution, and filter design for stationary and nonstationary 

random processes," IEEE Transactions on Signal Processing, vol. 58, no. 8, pp. 4079-4092, 2010. 
[13] F. Neeser and J. Massey, "Proper complex random processes with applications to information theory," IEEE Transactions 

on Information Theory, vol. 39, no. 4, pp. 1293-1302, July 1993. 
[14] B. Picinbono, "Second-order complex random vectors and normal distributions," IEEE Transactions on Signal Processing, 

vol. 44, no. 10, pp. 2637-2640, October 1996. 
[15] A. Papoulis and S. Pillai, Probability, random variables and stochastic processes, 4th ed. McGraw Hill Higher Education, 
2002. 

[16] S. Pei, M. Yeh, and C. Tseng, "Discrete fractional Fourier transform based on orthogonal projections," IEEE Transactions 

on Signal Processing, vol. 47, no. 5, pp. 1335-1348, May 1999. 
[17] J. Proakis, Digital communications. McGraw-hill, 1987, vol. 1221. 



DRAFT 



16 



[18] M.-H. Y. Soo-Chang Pei and T.-L. Luo, "Fractional Fourier series expansion for finite signals and dual extension to 
discrete-time fractional Fourier transform," IEEE Transactions on Signal Processing, vol. 47, no. 10, pp. 2883-2888, 
October 1999. 

[19] D. Tse and P. Viswanath, Fundamentals of wireless communication. Cambridge University Press, 2005. 
[20] B. Rimoldi, Course Notes: Principles Of Digital Communications. 

[21] H. Sun, G. Liu, H. Gu, and W. Su, "Application of the fractional Fourier transform to moving target detection in airborne 
sar," IEEE Transactions on Aerospace and Electronic Systems, vol. 38, no. 4, pp. 1416-1424, October 2002. 

[22] Y. Qin, L. Wenyao, Z. Shouli, and H. Hairong, "Detection of chirp signal by combination of kurtosis detection and filtering 
in fractional fourier domain," in 2nd International Congress on Image and Signal Processing CISP, October 2009, pp. 
1-6. 

[23] R. Jacob, T. Thomas, and A. Unnikrishnan, "Fractional Fourier transform based chirp detector versus some conventional 
detectors," in International Symposium on Ocean Electronics (SYMPOL). IEEE, November 2009, pp. 56-65. 

[24] I. S. A. Yetik and A. Nehorai, "Beamforming using the fractional Fourier transform," IEEE Transactions on Signal 
Processing, vol. 51, no. 6, pp. 1663-1668, June 2003. 



DRAFT 



